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The work is devoted to comparison of two different approaches to calculation of three- 
body resonances on the basis of the Faddeev differential equations. The first one is the 
well known complex scaling approach. The second method is based on an immediate 
calculation of the zeros of the scattering matrix continued to the physical sheet. 

For years the complex scaling [l],[J is considered as one of the most convenient methods 
for calculation of resonances. This method gives a possibility to rotate the continuous 
spectrum of an iV-body Hamiltonian with analytic potentials. The binding energies stay 
fixed in this rotation while the resonances appear as additional complex eigenvalues of 
the scaled Hamiltonian situated in the sectors of the unphysical sheets which become 
accessible during the scaling transformation. Surveys of the literature on the complex 
scaling method and its various applications can be found in ||. Here we only note that 
for a rather wide class of the interaction potentials the complex-scaling resonances are 
proved to be at the same time the scattering-matrix resonances Q (i. e. they are also the 
poles of the scattering matrix analytically continued in corresponding unphysical sheets). 

Various different methods, of course, are also used for calculations of the resonances. 
Among the methods developed to calculate directly the scatering-matrix resonances we, 
first, mention the approach based on the momentum space Faddeev integral equations f| 
(see, e.g., Ref. || and references cited therein). In this approach one numerically solves 
these equations continued into an unphysical sheet and, thus, the three-body resonances 
arise as the poles of the continued T-matrix. Another approach to calculation of the 
scattering-matrix resonances is based on the explicit representations |7|] for the T- and 
S-matrices, continued into unphysical sheets, in terms of the physical sheet. From these 
representations one infers that the three-body resonances can be found as zeros of certain 
truncations of the scattering matrix only taken in the physical sheet. Such an approach 
can be employed even in the coordinate representation 

To the best of our knowledge there are no published works applying the scaling approach 
to the Faddeev equations. Thus, the purpose of the present work is two-fold. On the one 
hand, we make, for the first time, the complex scaling of the Faddeev differential equations. 
On the other hand we compare the complex scaling method with the scattering-matrix 
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approach We do this making use of both the approaches to examine a model system 
of three bosons having the nucleon masses and the three-nucleon (nnp) system itself. 

We start with recalling that the three-body Schrodinger operator H 3b reads after the 
scaling transorm as follows 

= U{d)Hz b U{-d) = -e- 2 ' d A x + J2M^K\), a = 1, 2, 3. (1) 

a 

where d = i9 with S 6 I while x a ,y a are the standard Jacobi variables, X = (x a ,y a ). 
By Ax we understand the 6-dimensional Laplacian and by v a , the two-body potentials. 
The corresponding scaled Faddeev equations that we solve have the following form: 

[-e-% + v a (e*\x a \) - z]QW(X)+v a (e*\x a \)^QW(X) = f a (X), a = 1,2,3. (2) 

Here / = (f%, / 2 , ^3) is an arbitrary three-component vector with components f a belonging 
to the three-body Hilbert space. The partial-wave version of the equations @ for a system 
of three identical bosons with L = reads 

e-^H^^x, y) - z y) + V(e ie x)%(x, y) = f t (x, y), (3) 

Here, Hq denotes the partial kinetic energy operator and the partial-wave component 
of the sum m = $(i) + $(2) + $(3)_ Explicit expression for H{p as well clS db formula explicitly 
relating $^ to can be found, e.g., in ||. 

For compactly supported inhomogeneus terms fi(x,y) the partial- wave Faddeev com- 
ponent $i(x, y) satisfies the asymptotic condition 

y) = 8i ipd(e i6 x) exp{iy/E t - e d e w y) [a + o {y~ 1/2 )] 
exp(i x /E; e w p) r , _ 1/a v 1 

+ -j= [A{y/x) + o{p 1 )\ , 

where ipd{x) stands for the wave function of the two-body subsystem while ao and Ai(y/x) 
the main asymptotical coefficients effectively describing the contributions to $z from the 
elastic (2 — > 2) and breakup (2 — > 3) channels, respectively. 

In the scaling method a resonance is looked for as the energy z which produces a pole 
to the function 

<t>(6 1 z) = ([H F (6)-z}- 1 fj) (5) 

where Hp{9) is the non-selfadjoint operator resulting from the complex-scaling trans- 
formation of the Faddeev operator. This is just the operator constituted by the 1. h. s. 
parts of Eqs. (H). The resonance energies should not, of course, depend on the scaling 
parameter 9 and on the choice of the terms fi(x,y). 

In the table we present our results obtained for a complex-scaling resonance in the 
model three-body system which consists of identical bosons having the nucleon mass. To 
describe interacation between them we employ a Gauss-type potential with a barrier term 
taken from Ref. || the barrier value V& (see |§ for details) being fixed, V b = 1.5. The 
figures in the table correspond to the roots of the inverse function calculated 
on the basis of the finite-difference apporoximation of Eqs. (H|-fD in polar coordinates. 



(4) 
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Description of the finite-difference algorithm used can be found in Ref. |J. In the present 
calculation we have taken up to 400 knots in both hyperradius and hyperangle variables 
while for the cut-off hyperradius we take 40 fm. One observes from the table that the 
position of the resonance depends very weakly on the scaling parameter 9 which confirms 
a good numerical quality of our results. 



9 


^res (MeV) 


9 


Zres (MeV) 


0.25 


-5.9525 - 0.4034 i 


0.50 


-5.9526 - 0.4032 i 


0.30 


-5.9526 - 0.4033 i 


0.60 


-5.9526 - 0.4033 i 


0.40 


-5.9526 - 0.4032 i 


0.70 


-5.9526 - 0.4034 i 



We compare the resonance values of the table to the resonance value z res = —5.952 — 
0.403 i MeV obtained for the same three-boson system with exactly the same potentials 
but in the completely different scattering-matrix approach of Ref. ||. We see that, indeed, 
both the complex scaling and the scattering matrix approaches give the same result. 

As to the nnp system where we employed the MT I — III [|K| model, both these methods 
give no resonances on the two-body unphysical sheet (see ||). Also we have found no 
resonances in the part of the three-body sheet accessible via the scaling method. Thus, at 
least in the framework of the MT I — III model we can not confirm the experimatal result of 
Ref. |12| in which the point —1.5 ± 0.3 — i(0.3 ± 0.15) MeV was interpreted as a resonance 



corresponding to an exited state of the triton. 

The triton virtual state can be only calculated within a scattering-matrix method but 
not in the scaling approach. Our present improved scattering-matrix result for the triton 
virtual state obtained within the method |7||| is —2.690 MeV (i.e. the virtual level lies 
0.47 MeV below the two-body threshold). This result has been obtained with the MTI- 
III potential on a grid having 1000 knots in both hyperradial and hyperradial variables 
and with the value of cut-off hyperradius equal to 120 fm. Some values for the energy of 
virtual state obtained by different authors can be found in || , and all of them are about 
0.5MeV below the two-body threshold. 
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